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ABSTRACT 

Context. The observed distribution of orbital periods of Algols with a B-type primary at birth agrees fairly well with the prediction 
from conservative theory. Conservative evolution fails, however, to produce the rather large fraction of Algols observed with a high 
mass-ratio, especially: q e [0.4-0.6]. 

Aims. In order to keep Algols for a longer time with a higher mass-ratio without disturbing the distribution of orbital periods too 
much, interacting binaries have to lose a significant fraction of their total mass without losing much angular momentum before or 
during Algolism. We propose a mechanism that meets both requirements. 

Methods. In the case of direct impact the gainer spins up: sometimes up to critical velocity. Equatorial material on the gainer is 
therefore less bound. A similar statement applies to material located at the edge of an accretion disc. The incoming material moreover 
creates a hot spot in the area of impact. The sum of the rotational and radiative energy of hot spot material depends on the mass- 
transfer-rate. The sum of both energies overcomes the binding energy at a well defined critical value of the mass-transfer-rate. As 
long as the transfer-rate is smaller than this critical value RLOF happens conservatively. But as soon as the critical rate is exceeded 
the gainer will acquire no more than the critical value and RLOF runs into a liberal era. 

Results. Low-mass binaries never achieve mass-transfer-rates larger than the critical value. Intermediate-mass binaries evolve mainly 
conservatively but mass will be blown away from the system during the short era of rapid mass-transfer soon after the onset of RLOF. 
We have calculated the evolution of binaries with a 9 M© primary and a 5.4 M© companion over a range of initial orbital periods, 
covering case-A RLOF. Mass-loss from the system is achieved during direct impact only. 

Conclusions. We find systems that show Algolism for more than ten million years. RLOF occurs almost always conservatively. Only 
during some 20,000 years the gainer is not capable of grasping all the material that comes from the donor. The mass-ratio q e [0.4-0.6] 
which was hardly populated by conservative evolution now contains Algols for a significant fraction of their existence. 
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1. Introduction 

Eggleton ('2000') introduced the denomination liberal to distin- 
guish binary evolution with mass and subsequent angular mo- 
mentum loss from the conservative case, where no mass leaves 
the system . Liberal evolution must be at work since Refsdal et 
al. (119741) found no progenitor that can evolve into AS Eri in a 
conservative way. IMassevitch & Yungelson ( 1975 ) showed that 
in order to obtain agreement between theory and observations 
systems with combined mass < 6 M© have to lose 40 to 50 % of 
the mass lost by the donor. Sarna ( 1993 ) showed that only 60 % 
of the mass lost by the loser of JS Per was captured by the gainer, 
while 30 % of the angular momentum was lost during Roche 
Lobe Overflow (RLOF). It has been shown by the Brussels 
group that conservative calculations produce almost no Algols 
with large mass-ratios during case-B RLOF (Van Rensbergen, 
12003]) . Van Rensbergen et al. (l2QQ6l) included case-A RLOF 
into the comparison between observation and conservative bi- 
nary evolution theory, leading again to a too small number of 
Algols with large mass-ratios. The observed distribution of or- 
bital periods of Algols is, however, fairly well reproduced by 
conservative evolution. All the conservative evolutionary tracks 
can be found at http : //www . vub . ac . be/astrof ys/. In this 
paper we propose a scenario wherein violent phases of rapid 
RLOF can trigger mass-loss from the system. Whereas Algols 
at present go mainly through quiet phases of RLOF, they may 
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have had a violent past in which they have lost a considerable 
fraction of their mass. 



2. Observed orbital periods and mass-ratios 

To avoid an eternal confusion we replace the indices (1 = 
primary and 2= secondary) characterizing binary components 
by d and g. We use the index d for the donor, i.e. the star that 
will be the donor once RLOF has started, whereas g is used for 
its gaining companion. We define the mass-ratio q throughout 
as: 



De Loore & Van Rensbergen (I20Q5I) introduced the qualifi- 
cation Algolism for binaries during their Algol stage, when the 
mass of the donor has become necessarily smaller than the mass 
of the gainer. In our comparative study the value of the mass- 
ratio q will thus always be in the interval [0-1]. 

In this paper we compare theory with observations for 
Algol binaries with a B-type primary at birth. The catalogue of 
Budding et al. (2004 ) extended with the semi-detached Algols 
from Brancewicz et al. (1996 ) supplies us with 303 Algols which 
can be issued from conservative binary evolution with a B-type 
primary at birth. 
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The observed distribution of orbital periods and mass-ratios 
for the SB 2s among these systems is well established. The over- 
all distribution of mass-ratios also includes the SB Is and has 
been revisited since Van Rensbergen et al. ( 2006 ) claimed that 
more than 70% of the observed Algols are located inq e [0.4-1] 
if one uses qMs which is determined so as to make the param- 
eters of the most massive star fit main sequence characteristics. 
Two other methods evaluating mass-ratios of SB Is were, how- 
ever, not used in this study: the mass-ratio qic which is obtained 
by the light curve solution and qso, which uses the assumption 
of a semi-detached status. 

In this study we compared the values of qMs^ ^lc and 
qs D w ith the ^-values as determined by Pourbaix et al. 
(120041 http://sb9.astro.ulb.ac.be/). We found the ^5'/)-values not 
representative. The Pourbaix- values are underestimated by qic 
and overestimated by qMS • The observed mass-ratio distribution 
of the 303 Algols cited above has hence been recalculated using 
the mixing of qMs and qic that represents the numbers in the 
Pourbaix catalogue best, which leaves still 45% of the observed 
Algols inq e [0.4-1], as can be seen in Figure |2] 

3. Need for liberal evolution 

The conservative simulation has been explained in detail by Van 
Rensbergen et al. (2006). Figure [T] compares the observed or- 
bital periods of 303 Algols with a B-type primary at birth with 
the orbital periods obtained from conservative binary evolution. 
Larger initial periods leading to case-B RLOF produce mainly 
Algols with long orbital periods. Cases- A follow the observed 
distribution better. Since among the Algol population there are 
far more systems undergoing RLOF-A (fraction of the nuclear 
time-scale) than RLOF-B systems (fraction of the much shorter 
Kelvin-Helmholtz time-scale) the observed distribution of or- 
bital periods meets the results from conservative binary evolu- 
tion well. 

Figure [2] compares the observed mass-ratios of 303 
Algols with a B-type primary at birth with mass-ratios from 
conservative binary evolution. Cases-B produce more than 80% 
Algols with q e [0-0.2]. Cases-A produce most of their Algols 
with q G [0.2-0.4]. The fact that ^ 45 % of the observed Algol 
systems are in q e [0.4-1] forces us to state that conservative 
evolution can not be the only channel for the evolution of bina- 
ries with a B-type primary at birth. 

In our search for a liberal scenario we divide interacting bi- 
naries into direct impact systems and systems with an accretion 
disc around the gainer. During their evolution, binaries change 
positions in the (r - ^)-diagram, where q is the mass-ratio as de- 
fined by relation ([T]) and r the relative radius of the gainer: i.e. its 
radius divided by the semi major axis of the system. Using the 
semi-analytical ballistic calculations of Lubow and Shu (1975 ), 
two curves co^ and ojmin are drawn onto the (r - ^)-diagram. If a 
system is located above ojd the gas flow coming from the donor 
will hit the gainer directly. If the relative radius of the gainer is 
below comin, the gas stream will feed a classical accretion disc 
with relative outer radius equal to 2 x 6i;j. We performed our 
calculations so that systems between oj^ and ojmin will develop 
an accretion disc with a relative outer radius growing gradually 
from ojd = Rg to 2 X oj^. 

Furthermore, the degree of liberalism will be measured with 
the quantity JS e [1-0] which is defined as follows: 

Mg =-l3Mf^^ (2) 



^RLOF • 



is the negative value of the mass lost by the donor and 
Mg the positive value gained by the gainer, (l-p) M^^^^ is hence 



the negative value blown away by the gainer. J3 = I characterizes 
conservative evolution, where every amount of mass lost by the 
donor is captured by the gainer. Figure [6] shows an evolution 
which is conservative most of the time, but in which up to ^ 80 
% (i.e. p ^ 0.2) of the mass lost by the gainer through RLOF is 
blown into interstellar space during a liberal era. 

4. Liberal evolution during direct impact 

4.1. The rapidly rotating gainer 

Conservation of angular momentum spins the gainer up due to 
the impact of RLOF-material coming from the donor star. Mass 
located near the gainer's equator gets loosely bound when the 
gainer rotates rapidly. The spinning-up of the gainer is character- 
ized by an enhancement of its rotational angular momentum AJ^ 
which is given in cgs-units by Packet ( 1981 ), corrected with the 
impact-parameter d as shown in figure |4] simulating the direct- 
hit scenario (Langer J1998l) : 



M, AM, 



AJ^ = 6.05 X 10^^ X [^]2 X [^ + 



AM, 
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Mc 



2Mn 
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X — (3) 

Mo Rg^ ' 



This spin-up is, however, counteracted by tidal interactions 
which were first studied by Darwin ( 1879J. The formalism for 
tidal downspinning can be talcen from Zahn (I1977K who gives a 
suitable approximation for the synchronisation time- scale: 



Tsync (jr) = q X [ — ] 
^8 



(4) 



This expression uses the semi major axis a of the binary and 
a mass-ratio q, in which the star that has to be synchronized is in 
the denominator. This is the gainer in our case, so that ^ = ^ • 

Tidal interactions modulate the angular velocity of the gainer 
ojg with the angular velocity cogrb of the system. According to 
Tassoul (2000) one can write: 



1 



OJg - CO orb 



dCDorb 

dt 
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' sync 



sync 
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''sync 



(5) 



Using the moment of inertia Ig of the gainer we find the ex- 
pression which was used by Detmers et al. ( 2008 ) in his scenario 
for the liberal evolution of a massive close binary: 



AJ = IgX (Ojorb - ^^^) X [1 - ^ 



' sync >^ J sync 



(6) 



Tidal interactions spin the gainer down when ojg > coorb- 
Tides spin the gainer up when cOg < oJorb- fsync - 1 represents 
weak tidal interactions whereas fsync = 0.1 implies strong tides. 

When the upspinning stops at the end of RLOF, tidal in- 
teractions will settle the system into a situation with cOg = 
(jJorb- Expression ([6]) then implies that AJ~ = AJ^ = 0. 
Synchronisation is achieved and angular momentum remains 
conserved. 

The spin-up of the gainer is strongly at work during the 
short era of rapid mass-transfer soon after the onset of RLOF. 
Disregarding tidal interactions. Packet ( 1981 ) showed that crit- 
ical velocity at the gainer's equator is achieved more easily if 
the mass of the gainer is small. Figure [3] shows the evolution 
with time during the rapid phase of mass-transfer of a (9+5.4) 
Mq binary with an initial orbital period of 3.6 d. The spin-up 
works in this case during hydrogen core burning of the donor. 
From Figure [3] it is clear that critical velocity is only achieved 
assuming weak tidal interactions, whereas strong tides prevent 
the gainer to rotate at critical velocity. 
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Fig. 1. Observed distribution of periods of Algols compared to conservative evolution of binaries with a B-type primary at birth. 
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Fig. 2. Observed distribution of mass-ratios of Algols compared to conservative evolution of binaries with a B-type primary at birth. 
The fact that ^ 45 % Algol binaries are observed in ^ € [0.4-1] shows that conservative evolution can not always be valid. 
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Fig. 3. The gainer in a (9-^5.4) Mq binary with an initial period of 3.6 days rotates at critical velocity during the rapid phase of 
mass-transfer if the tidal friction is weak. Strong tidal interaction prevents very fast rotation of the gainer. 
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4.2. The hot spot 

4.2.1 . Visibility of a hot spot 

The accretion luminosity Lace on the gainer's surface is caused 
by the impact of RLOF-material. This impact causes a bright 
spot a t the trailing side of the gainer's equator (Peters & Polidan, 
I2QQ41) . Such a hot spot can also be created at the outer edge of 
an accretion disc, leaving the possibility for material to spiral in- 
ward. The geometry of the direct impact system shown in figure 
|4 illustrates that the hot spot is only visible to the observer near 
phase O ^ 0.75. This is the case as well if the hot spot is located 
at the gainer's equator or at the edge of an accretion disc. 

4.2.2. Characteristics of a hot spot 

When matter falls from infinity into the potential well of the 
gainer it gets hot on the gainer's surface and emits radiation. 
The accretion luminosity is then given by the diff'erence of the 
potential energy at infinity and at the point of impact P. 



LZc = U(oo) - U(P) = Gx 



(7) 



This accretion luminosity is weakened because matter hits 
the gainer coming from the first Lagrangian point Li which is 
diff'erent from the idea that it starts from infinity. Taking into ac- 
count the sonic speed in Li, a reduction factor D e]0 - 1[ can 
easily be calculated using the appropriate distances in the coro- 
tating system shown in Figure |4] The real accretion luminosity 
can then be written as Lace = D x and can be calculated 
numerically as: 



^ = U(Li) - U(P) = 3.14 X 10^ X Z) X 

Lq 



Mo 



(8) 



The luminosity of a star is a global quantity. The accretion 
luminosity Lace, however, has to be evaluated locally. If we want 
matter to escape from the gainer's equator, it has to be removed 
from the restricted surface area of the accretion zone, which is 
smaller than the entire gainer's surface. Hence Lace has to be 
evaluated as a local quantity which is strengthened by the lim- 
ited area of the accretion zone in which the accretion luminosity 
is concentrated: ^ can be used as the outward radiative pres- 
sure in the area of impact. S is the fraction of the stellar surface 
covered by the accretion zone (S <^ 1). 

On the other hand. Lace is weakened by the low efficiency of the 
accretion luminosity and only ^ (A » 1) can be converted into 
radiative flux. 

The radiation pressure exercised by a hot spot is thus pro- 
duced by : 



S xA 



K 



with K = S xA 



(9) 



The numbers A and S (and hence the crucial quantity K) can 
be derived from the observations using the following procedure. 

When the gainer has no spots its luminosity is L^. The luminosity 
of the spotted gainer is then given by: 



L^ T 

Lq Lq 



(10) 



Ladd is the fraction of the accretion luminosity Lace, given by 
relation (ISl), which is converted into radiative flux. 



• Ladd < implies a dark spot 

• Ladd = implies a non spotted stellar surface 

• Ladd > implies a hot spot. 

We assume that the hot spot is created by RLOF-material 
infalling from the donor star. 

Next, we can rewrite equation ([TO]) as: 



8 _ p ^g ^spot ^ eff,g _^ ^ spot ^ ^ spot 

Lq~^ Sq (5770)4^ ^ (5770)4^ 

For a hot spot with Tgpot > Teff,g we can write this as: 



'-^add ^ spot 



1 



X 



Lq Sq (SllOy 



X l-T^pot ^tffj 



(11) 



(12) 



Knowledge of Ladd (e.g. through direct observation) deter- 
mines the quantity A as: Ladd - Since Ladd ^ ^aee, we have 
that A » 1. 

Temperature and size of the spot are related through relation 
([12]), which can be transformed into: 



spot 



+ X ^ X (5770)4] 

^ spot 



(13) 



Every measured value Lace is reproduced by an infinite num- 
ber of combinations of S spot and Tspot- Only a restricted range 
of values hereby reproduce a realistic hot spot. Knowledge of 
(e.g. through direct observation) determines the quantity 



^ spot 

S as 



Since S , 



<^ S we have that S «: 1. The crucial 



Sg • ^^^^^^ ^ spot 

quantity K, as defined by equation ([9]), is thus the product of a 
factor A » 1 and a factor 5" ^ 1 . It is clear that small values of K 
will support mass-loss whereas large values of K will suppress 
mass-loss from the system. 

It has to be noticed that the quantity K can be evaluated di- 
rectly from equation ([12]) which can be written as: 



KxLq 



Rg ' 1 



] x[r: 



spot 



' eff,g^ 



(14) 



One can thus consider /sT as a function of Lace and Tspot, a 
relation that determines the important quantity K without dis- 
entangling it into its constituents A and S . Measurement of the 
temperature of the spot and determination of the mass-transfer- 
rate (often using observed changes of the orbital period), which 
determines the accretion luminosity through relation ([S]) enable 
us to determine the quantity K throughout the entire evolutionary 
computation. 



4.3. Mass can leave the system 

If one thinks about possible mass-loss from the gainer into the 
interstellar medium we first have to evaluate the binding energy 
of a test mass m in the hot spot point P of Figure]?] Without ro- 
tation of the gainer and without accretion luminosity the energy 
of a test mass m at rest in P is given by: 



U = - G [ + J - - mr [ J 



Sd 



(15) 



orb 



This negative amount of energy is counteracted by the lu- 
minosity of the gainer (Lnuei), its enhanced equatorial velocity 
and the positive energy input Lace given by relation ([5]). The im- 
pacting gas stream has, however, already been used to spin up 
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Fig. 4. Geometry determining the binding energy of a test mass in the impact point P in the middle of a hot spot. The observer looks 
into the direction of the small arrows. The hot spot P is turned towards the observer only near phase O ^ 0.75. 



the gainer as explained in section ( |4.1| ). The rotational kinetic 
energy is hereby raised at a rate Krot- The energy rate which re- 
mains available to build the hot spot is thus a little less than Lace 
and equal to: 



^nrr ^aCC 



In this paper we have replaced k by kross- The Rosseland 
opacities at the gainer's surface for solar abundances and the ap- 
propriate values of log T^ff and log g were hereby taken from 
Kurucz (11979) , defining the quantity F as: 



(16) F 



KRo 



(20) 



This remaining part {L~^c) ^^^^ 1^ increase the temperature 
of the hot spot. When the rotational energy ^ m v^^ of a test mass 
located at the gainer's equator approaches the break-up velocity 
already almost 50 % of the binding energy given by relation ( 15 ) 
is surmounted. 

The quantity K rot is n ot easy to observe and has thus been ne- 
glected in section (4.2.2) but is easily followed up during a run 
of the Brussels binary evolution code. When the gainer rotates 
rapidly, the spinning up becomes almost impossible and practi- 
cally all the energy available from the accretion will then be used 
for building the hot spot, so that L~^^ ^ Lace- 

The work done by the outward radiative force is: 



rad 



m K Lnucl + 

~c~ 47iR 



iKcc) 



(17) 



One obtains the total amount of energy on a test particle at 
the gainer's equator by adding the rotational kinetic energy and 
the energy input as given by relation (17) as positive terms to 
equation (15). The test mass can only leave the system, when 
this energy is positive. So matter in the hot spot can only leave 
the system if (after convenient division by m): 



1 2r 27r 

2 Port 



j^RLOF,crit ^j^^^ ^^j^^ M^^^^ which causes a value of 



Lacccrit such that the left hand side of equation ( 18 ) equals zero. 
Every amount of M^^^^ in absolute value exceeding m^^^^^^^^^ 
will leave the system. In other words: supercritical mass will be 
lost by the system. 

For simplicity one could take only the Thomson scattering 
into account to evaluate relation ([18]): 



K = Ke = 0.34 cm g 



(19) 



Figure [5] shows that at low temperatures the Rosseland opac- 
ity is smaller than given by relation (19), whereas Thomson 
scattering produces the total opacity at high temperatures. 
From 11,000 K on, however, there is a temperature range 
wherein the Rosseland opacity is the largest, so that it de- 
mands less energy for the test particle in the hot spot to 
leave the system. Since the temperature of the hot spot ex- 
ceeds the temperature of its surroundings, we calculate in our 
Brussels binary evolutionary code the quantity F at a tempera- 
ture equal to 1.5 x Te/f^g. 

With these assumptions, a numerical value for m^^^^^^^^^ can 



be found from equations ( 18 ) and m: 



I^^LOf.c„Y| ^ ji 224 lO-'[— + — ] 

' ' FDMg^ Sd Rg 

+ 8.219 10"^ l—f - 3.208 10"^ 



orb 



+ 3.183 10" 



Kr, 



3.183 10-« — L„„rf} 



(21) 



Every quantity is expressed in the usual units and Veq in 



km 



Mass-loss out of a binary is mainly driven by rapid rotation 
and large mass-transfer-rates for systems with a B-type primary 
at birth. For very massive binaries, however, Lnud may be suffi- 
ciently large so as to drive matter out of a binary in which the 
gainer is spinning rapidly (Petrovic et al. .2005j . 

5. Possible liberal evolution with an accretion disc 

The geometry of a binary with an accretion disc was introduced 
in section ([3|, yielding accretion discs with outer radius Rdisc be- 
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Fig. 5. Opacity ratio as a function of temperature, showing the dominance of kross in a temperature region above 1 1,000 K. 



tween Rg = ojd x a and 2 x aj^ x a. Replacing Rg with R^isc 
relations ^ and ([S]) can be used to calculate the accretion lumi- 
nosity in the point P at the outer edge of the spot. This accretion 
luminosity is considerably reduced when the distance between 
Li and P is small. 

For the calculation of the critical amount of mass-transfer. 



forcing the system to evolve in a liberal way, equation ( [21] ) has to 
be modified. Dealing with a Keplerian accretion disc, the equa- 
torial velocity of the gainer has to be replaced by the Keplerian 
velocity of a test mass located at the edge of the disc. Due to 
the virial theorem, the rotational kinetic energy of the test mass 
covers exactly 50 % of the gravitational potential energy caused 
by the gainer. In this way, one can approximate the quantity 

j^RLOF,crit ^.^j^. 



^^RLOF,crit^ 



KRl 

_ disc J 

F DMg 
+ 8.219 10" 



[1.224 10-^[— + 



2xR, 



disc 



orb 



(22) 



Our calculations as summarized in section ^ show that ini- 
tially (9+5.4) Mq binaries evolve almost always conservatively, 
except during a short era of rapid mass-transfer short after the 
onset of RLOF. This situation occurs always during hydrogen 
core burning of the donor and with the direct impact geometry. 
The accretion disc is only formed when the orbit widens. The 
mass-transfer-rate is, however, not sufficiently large in this case 
so that the value given by expression ([22]) is never attained. With 
our simple assumptions, disc systems that are formed after an era 
of RLOF-A, evolve conservatively. Intermediate-mass binaries 
with larger initial orbital periods so that RLOF-A is avoided and 
an accretion disc is formed (not included in this paper) might 
very well also experience an era of liberal evolution. 

The mass-loss rates in Table 1 cover the direct impact 
systems only. Mass-loss in the case of an accretion disc 
around the gainer would enhance the degree of liberalism men- 
tioned in Table 1. Subsection ( [7.2| ) mentions, however, sev- 
eral intermediate-mass binaries with an accretion disc by which 
mass-loss from the system was observed. Hydrodynamical cal- 
culations by Bisikalo et al. (120051 ) follow matter that spirals in- 
wards and creates a subsequent hot line, caused by the shock 
issued from the interaction of the circumstellar disc with the 



flow from the inner Lagrange point Li . This concept was intro- 
duced in three-dimensional numerical simulations by Sytov et al. 
(120071) for a system with an orbital period of 0.23 d, a low-mass 
donor of 0.56 Mq and a 0.6 M© white dwarf gainer. Even this 
low-mass system with an accretion disc around the gainer loses 
mass at a rate of 3 x 10"^^ ^ through the third Lagrangian 
point L3 . This type of mass-loss is not yet included in our com- 
putational code. 



6. Calibrating the quantities A and S 

In this section we will evaluate relation ( [2T] ) carefully. 

• D = ^ is followed up during the entire evolutionary 

^acc 

calculation 

• F = is given introducing Tspot into figure |5] 

The quantity K is, however, more difficult to determine. For 
only a few binaries observations of Lace, Ladd, S spot are available 
in order to determine the quantities A and S and subsequently K: 



Lace 
Ladd 



• s = ^ 



K = AxS 



The situation where Tspot has been observationally deter- 
mined is largely to be preferred, because the quantity K can then 
be directly evaluated using relation ( [T4] ), which can also be writ- 
ten as: 



K 



X (5770)4 



(23) 



It is clear that the value of K using relation ( 23 ) can diff'er 
from the {K = Ax 5" )-calculation because all observed stellar 
parameters are not always entirely internally consistent. 
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6.1. Gauging the quantity S 

Gunn at al. J1999I) use an expression of Pringle (11985b giving 
the area of the cross section of the stream starting from the first 
Lagrangian point Li with the local sonic speed v sonic towards the 
gainer, for the calculation of the surface area of a hot spot around 
the impact point P: 



The conservative relation ( [25] ) has been used in section ^ 
because most of these systems are in a quiet phase of conserva- 
tive mass-transfer. The strict application of the conservative limit 
in liberal cases underestimates the real amount of mass-transfer 
only very slightly (Erdem et al. 2007 ), the diff'erence being well 
below the uncertainty introduced by the (O-C)-procedure. 



' spot 



2.77 10^^ T~' 



\LiP\ R^3 



(24) 



In this expression (24) we have used v sonic = 0.125 ^/f 
The number 0.125 has oeen taken between 0.1 and 0.15, num- 
bers which are valid using solar abundances for respectively 
a neutral gas and a completely ionized plasma. The quantity 
S = can now be calculated for every binary. Surface ar- 
eas of the spot have been determined for a few binaries: VW 
Cep ( Pustylnik & Niarchos 2000), CN And (Van Hamme et al. 
l2QnT1) , KZ Pav (Budding et al. 2001 ) and V505 Sgr (Gunn et al. 
I1999K showing a fair agreement with the results from expression 
(|24l). 



7. Calibration of tlie quantity K 

In this section we list the few semi-detached binaries for which 
one or more of the quantities A, S and K have been determined 
within reasonable limits. Since K = A x S, two of the three 
quantities allow to determine the third one. Simultaneous mea- 
surement of the three quantities is, however, very useful in or- 
der to narrow the error bars around the individual estimates. We 
remind the reader that the quantity A is known from Ladd^ the 
quantity S from S spot and the quantity K from Tspot- It is obvious 
that more accurate and a greater number of observations would 
highly improve the statistics outlined below. Especially accurate 
determinations of spot temperatures would make a more precise 
determination of the crucial quantity K possible. We have sorted 
the systems by ascending values of their total mass. 



6.2. Gauging the quantity A 

Section ^ lists the semi-detached binaries for which ob- 
servations are available allowing to determine some or 
all of the quantities A, S and K. The cases for which 
S(S spot)^ MLadd) (^nd K(Tspot) were measured or determined 
separately but for which the relation K = Ax S was obviously 
violated are not included. 

The quantity A is defined from the measured quantity 
Ladd and the calculated quantity Lace, as A = Lace 
is determined in our sample of semi-detached binaries in 
section ([Tj) using the values of ^ which were determined 
with the widely used (0-C)- procedure explained in detail 
by Sterken (120051) . Observed times of light minima (O) of 
eclipsing binaries were taken from the website of Kundera 
(120021 http : //www . oa . u j . edu . pl/ktt/kr ttk_dn . html), 
completed for the years after 2002, with the data published in 
various issues of the "Information Bulletin on Variable Stars, 
Commission 27 of the I.A.U.", which we found for each binary 
at http://simbad.u-strasbg.fr/simbad/. Subsequently, a 
best fit parabola was drawn through the observed times of light 
minima (C). The (0-C)- procedure then yields the observed 
value of ^ (expressed in in ^). The parabolic fit is sometimes 
convincing and sometimes doubtful. The latter cases have been 
disregarded because they occur when mechanisms diff'erent from 
mass-transfer (e.g. magnetic braking) are at work. 

When a convincing parabolic behavior ^ is available, the 



negative value of ^ (expressed in in — ) was calculated with 



Me 



the conservative relation J25|. We use this value of ^ to cal- 
culate Lace with relation 



dP 
~P~ 



^dMd ^dM, 

m7 



+ 3- 



(25) 



This equation is the conservative limit of the more general 
relation including mass-loss (dM 9^ 0) and angular momentum 
loss (dJ 9^ 0) from the system. 



dP ^dMd dM, 
— = 3 +3- 



Md 



_^dJ _ dM 



(26) 



7.1. Direct impact systems 
7.1.1. VW Cep 

Stellar parameters and restricted ranges for Ladd, S spot and Tspot 
from Pustylnik & Niarchos ( 2000 ), internally consistent values 
of f and ^ from Devita et al. ( 1997 ) and Pribulla et al. (»2000l) 
determine the value of K in the range [0.1-0.2]. 



7.1.2. CN And 

Stellar parameters from Van Hamme et al. (I2001L ^ from 
Samec et al. (1998), spot temperatures and sizes from Van 
Hamme et al. ( 2001 ) and references therein determine the value 
of K in the range [0.08-0.18]. 



7.1.3. KZ Pav 



Stellar parameters and approximate values of Ladd and ^ and 
thus ^ are given by Budding et al. ( 2001 ). Consistent values of 
^ and ^ can be found in Walker & Budding (1996 ). The value 
of S spot can been calculated with the expression ( 24 ) or can be 
calculated using the isomorphy with the very similar systems CN 
And and VW Cep for which S spot has been measured. Values of 
Tspot found from Ladd and S spot are comparable to those of CN 
And and VW Cep. We find a value of K in the range [0.03-0.08]. 



7.1.4. V361 Lyr 

Stellar parameters are given by Hilditch (12001b and Yakut & 
Eggleton (120051) . The mass-transfer-rate has been established by 
Hilditch (I2OQTI) as 2. 1 8 X 10"^ ^ , who also gives a value of Tspot 
= 9,500 K and ^ = 0.5. These numbers generate A = 2.6 and 
^=0.36. From this one finds S =0.14, which is more than one or- 
der of magnitude larger than the surface area of the spot given by 
relation pA\\ S = 0.011. Since relation ([24|) yields surface areas 



that often overestimate the real spot surface area, a spot surface 
area covering 14 % of the whole stellar surface is extremely high. 
With S = 0.011 we find K = 0.028 and Tspot = 17,700 K. Hence 
we have taken a mean value of K for calibration. 
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7.1.5. RTScI 



Stellar parameters and Tspot are given by Banks et al. (119901) . 
A narrow range of possible values of ^ and hence ^ can be 
found in Clausen & Gron bech O gTTj), Rafert & Wilson ( 19841) 
and Duerbeck & Karimie (119791) . Using these data we find K in 
the range [0.02-0.04]. With the value of S given by expression 
([24|) we also determine the efl&ciency factor A in the range [3-6]. 



7.1.6. U Cep 



Stellar parameters are taken from Budding et al. (|2001|) . 
Pustylnik (I1995t quotes a mass-transfer-rate of ^ 10"^ ^ but 

we will use the slightly smaller value 5 x 10"^ ^ as derived 
with the (0-C) method as outlined in section \62\ . A hot spot is 
reported from observations with i\\tFUSE spacecraft by Peters 
(120071) . The spot temperature is found to be ^ 30,000 K, lead- 
ing to ^ 3.3 X 10"^ for this system. The impact on the surface 
of the gainer creates a splash zone in which large velocities are 
identified in the spectrum. In order to keep the quantity A within 
physically possible limits, the impact zone should in this case 
be an order of magnitude smaller than S 0.02, as reported by 
Peters C2007i) . 



7.1.7. USge 

A consistent set of stellar parameters is given by Kempner & 
Richar ds (fT999l) . Richards & Albright ( 11^991) and Vesper et al. 
(120011) . Manzoori & Gozaliazl (i2007i) have determined a mass- 
transfer-rate of 1.79 X 10"^ A hot spot with a temperature 
in the range [20,000 - 100,000] K has been seen by Richards 
& Albright (1999). Hot circumstellar gas at temperatures in the 
range [60,000-200,000] K has been identified by Kempner & 
Richards (119991) . With a hot spot surface area given by relation 
( [24| ) only the lowest possible spot temperature of ^ 20,000 K can 
be accepted. Higher spot temperatures violate the law of conser- 
vation of energy, because they generate values of the quantity A 
below 1. Consequently we find a value of ^ 0.02. 



7.1.8. SV Cen 



Stellar parameters as given by Brancewicz & Dworak (119961) 
predict an orbital period increase rather than the observed de- 
crease. Stellar parameters from Wilson & Star (1 976) and 
Drechsel et al. (119821) yield contact systems in which there is 
no room for the construction of a hot spot. A small sized hot 
spot with a high temperat ure (T ^ 10^) has, however, been ob- 
served by Drechsel et al. (119821) in the UV. In this section we 
have used the stellar parameters of Rucinski et al. (,1992) which 
yield a close semi-detached binary. Herczeg & Drechsel (119851) 
notice that ^ is very unsteady giving a value of ^ in the range 



[1-4] X 10 



dt 

-4 Mq 



But even such high mass-transfer-rates cannot 



create a temperature of 10^ K in a hot spot with a size as given by 
relation (24). Depending upon the mass-transfer-rate we obtain 
spot temperatures in the range [30,000-45,000] K. The largest 
value of K that can be obtained is this way and will be used in 
this analysis is ^ 4.64 x 10"^. The spot temperature of ^ 10^ K 
would originate in a very small splash zone characterized by K 
^ 10"^ which according to relation (21 ) would lower the critical 
amount of mass-transfer below any reasonable value. 



7.2. Systems with an accretion disc 

In this case the hot spot is created at the edge of an accretion 
disc. The surface area of the spot is now hardly defined by re- 
lation (24). The quantity K can now be derived by comparing 
the temperature of the hot spot with the temperature of the edge 
of the disc rather than with the surface of the gainer. The radius 
of the accretion disc has been calculated as explained in section 
([3]). For the temperature profile of the disc we take the formulas 
as derived from first principles by Carroll & Ostlie (.2007 J : 



T(r) = T, 



disc 



' disc 




478, 074 [ 



(27) 



Masses and radii are in solar units and mass-transfer-rates in 
solar masses per year. The systems are again sorted by ascending 
values of their total mass. 



7.2.1. KO Aql 



Stellar parameters are taken from Soydugan et al. (120071) and 
Vesper et al. (2001). The data of Panchatsaram & Abhyankar 
(1981) yield a detached system and have therefore been disre- 
garded. ^ is taken from Panchatsaram & Abhyankar (11981! ) and 
Soydugan et al. (2007). This system has an accretion disc with 
Rdisc ~ 1.5 xRg. Expression (27 ) calculated with Mg and Rg 
yields a disc temperature of ^^,500 K which drops to ^ 4,600 
K at the edge. Soydugan et al. (2007) have observed hot spots 
in this system in which the gainer spins 30 % faster than syn- 
chronously. The quantity K can unfortunately not be estimated 
without measurements of Tspot- 



7.2.2. SW Cyg 

This system with an accretion disc shows much similarity with 
TT Hya. A consistent set of stellar parameters has been given by 
Richards & Albright (1999), Budding et al. ( 20041) . Brancewicz 
et al. (1996) and Vesper et al. (2001). Qian et al. (l2QQ2l) have 
determined = 2.15 x 10"^ — from a measured value of 

d t y 

^. Relation (pTj) yields Tdisc = 6,300 K and Tedge = 3,100 K. 
Albright & RicMrds (1996) found that the outer edge of the disc 
of TT Hya emits 5 times as much radiation than the same 
structure of SW Cyg. Gauging this with the spot temperature of 
TT Hya we find for SW Cyg: Tspot ~ 11,300 K. The value of K 
is in this case located around 0.016. 



7.2.3. TT Hya 

A coherent set of stellar parameters for this system is given by 
Eaton & Henry (1992), Peters & PoHdan (2004), Van Hamme 
& Wilson (1993), Richards & Albright (1999), Kulkarny & 
Abhyankar (1980) and Panchatsaram & Abhyankar (1981). The 
geometry of the system is thus well defined and a detailed hy- 
drodynamical model is given by Miller et al. (2007). Peters & 
Polidan ( 2004) quote a small value of Ladd and a spot temper- 
ature of 17000 K. Unfortunately there are no measures avail- 
able of ^ which makes a reliable determination of ^ diffi- 
cult. Peters & PoHdan (fT998l) argue a lower limit of 10"^^ ^, 
whereas Miller et al. (120071) determine a somewhat larger mass- 
transfer-rate of 2 X 10"^^ which agrees with the low value 
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of Ladd but contradicts the high spot temperature and t he dis c 
temperature which is evaluated at 7,000 K by Miller et al. (2007 ). 

A mass-transfer-rate of 2 x 10"^^ — would create an accretion 

y 

disc around the gainer of TT Hya of 1,350 K, leaving only 400 
K at its edge. The observed spot temperature of 17000 K would 
only be achieved by an extremely small sized spot character- 
ized by an impossible value of K 10" '^). If at the other hand 
the mass-transfer-rate would be as high as 2 x 10"^ ^ (as 
is the case for the similar system SW Cyg) one would obtain a 
disc temperature of ^ 7,000 K leaving still 2,500 K at its edge. 
The value of K would then have been reduced to a somewhat 
more acceptable value of ^ 10""^. This value will, however, not 
be taken into account for our calibration because the data on ^ 
(the quantity A) and Tspot (the quantity K) seem to contradict 
one another. 



7.2.4. V356 Sgr 

Stellar parameters are from Polidan (11989b , Simon (119991) and 
Peters & Polidan ((2004j). From the observed time series of ^ 

Polidan (J_989 ) determines a mass-transfer-rate of 4 x 10"^ 
The gainer rotates at 6.3 times the synchronous velocity (Simon, 
1199 9 ). Ladd has been measured by Peters & Polidan ( 2004 ) who 
also notice that a significant fraction of the circumstellar material 
is located near the surface of the gainer. This binary is in its era of 
mass-loss from the system so that the mass-transfer-rate exceeds 
its critical value given by relation (22), leading to a value of the 
quantity K of at least 1.9 x 10""^. 



include the quantity K in the Brussels binary evolutionary code 
we have used the following tentative best fit relation: 



K = 0.228 X [Md + Mg] 



-1.735 



(28) 



This relation has to be considered as very provisional. A sim- 
ilar but more accurate relation will, however, only be obtained if 
more precise and a greater number of observations on L^^j, S spot 
and especially Tspot will be available. 



8. Results 

A binary will lose mass only when the mass-transfer-rate rises 
above the critical rate as given by relation (21 ) for a direct im- 
pact system or relation ( [22| ) if the the hot spot is located at the 
edge of an accretion disc. A large mass-transfer-rate will thus be 
the major driving mechanism of mass-loss from the system. A 
sufficiently large rate is achieved during a short era of fast mass- 
transfer soon after the onset of RLOF. Such high mass -transfer- 
rates are predicted by the binary evolution codes of: 

• Paczynski 

• Eggleton 



Brussels which is used in this analysis 



Paczynski (11967b quotes a maximum value of ^ as high as 
3.4 X 10"^ ^ for a (16+10.67) M© binary with an initial period 
of 5.55 days. This value is almost recovered by the Brussels 
code in its conservative mode (1.8 x 10"^ ^). According to 
Paczynski (1967) this system transfers 9.1 M© from donor 
to gainer in 4,000 years only, during hydrogen shell burning of 
■ ^ ^ ^ ^ the donor. This amount of mass will certainly be too large for 

, , the gainer to be captured completely. Using the Eggleton code 

Harmanec & Scholz (11993]) determine values of ^ and ^ ([19711 & [20021) . Yungelson (200l finds a mass-transfer-rate up 
which are very high (^ = 3.4 x 10"^) ^ and are confirmed to 5 x 10"^ ^ soon after the onset of RLOF during the hydro- 
by Ak et al. (2 007 ). T his value is sHghtly larger than the one gen core burning of the donor for a (9+5.4) M© system with an 
used by Simon (I1999D. Lyr is a binary with 7^/^,^=28,000 K, initial period of 3 days. Also this mass-transfer-rate is almost 
Tdisc=^ fiOO K and 7,^^^=20,000 K (Harmanec, 2002). The best predicted during the era of fast mass-transfer in this system by 

the Brussels code in its conservative mode (2 x 10""^ ^). 

Low-mass binaries hardly achieve the requirements for sig- 
nificant mass-loss from the system. In this section we report the 
evolution of an initial binary of (9+5.4) M© with initial periods 
ranging from 1.4 to 5 days. We have calculated the evolution us- 
ing the weak as well as the strong tidal interaction as outlined in 
section [3] The quantity jS, defined by relation ([2]), was frequently 
treated as a free parameter. In this study we calculate jS through- 
out the entire evolution of the binary. 

Figure |6] shows the behavior of JS as indicator for the liberal 
evolution of a (9+5.4) M© binary with an initial period 3.2 
days. The two diff'erent tidal assumptions make hardly any dif- 
ference for the result. The binary will almost always evolve 
conservatively, but in a short and violent era soon after RLOF 
ignition, mass will be lost from the system: 2.54 M© adapting 
the strong tidal interaction and 2.74 M© when the tidal interac- 
tion is weak. The small diff'erence between the two results can 
be understood from relation ( [2T] ). The critical mass-transfer-rate 
is indeed lowered by a high equatorial velocity of the gainer, 
which is more easily achieved with weak tides. However, the 
same critical rate is increased by higher rates of rotational ki- 
netic energy of the gainer, which is favored by weak tidal in- 
teraction. In many cases both eff'ects more or less cancel each 
other. Notice in Figure |6] that the system frequently undergoes 



spot 

model is a thick Keplerian accretion disc around the gainer in 
which the eff'ective gravitation at the edge is 0. Evaluated at 
the edge of the accretion disc we find Tedge = 7,000 K, the accre- 
tion luminosity (428 L©) and K 0.012) evaluated with relation 
( [23] ). Harmanec (2002) notices the existence of bipolar jets with 
wind velocities comparable to wind velocities of O stars. /3 Lyr 
is blowing mass into the interstellar medium because its mass- 
transfer-rate is above the limit given by relation ([22]) which is 
valid for a hot spot on the edge of an accretion disc. In that case 
the quantity K needs to be ^ 0.0066 so that (^ = 3.4x10"^ ^) 



equals the critical value given by relation (22). For the calibra- 
tion of the quantity K we have taken the last value which would 
also have been obtained with Tspot = 23,000 K a number that 
diff'ers not significantly from the value of 20,000 as K quoted by 
Harmanec ([20021) . 



7.3. K as a function of total mass 

In this section we have determined values of K for eight systems 
undergoing direct impact and three systems with accretion discs. 
No particular trend is found for the quantity ^ as a function of 
many stellar parameters. However, one finds that low-mass bina- 
ries have larger values of K than more massive ones. In order to 
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two separate stages of mass-loss {fi <^ 1). The violent epoch 
of mass-loss indeed starts with Md > Mg. The orbital period 
will shrink until = Mg, which makes the impact parameter 
d (Figure |4]) so small that the accretion luminosity as calculated 
with expression ^ will drop below the critical value. As soon 
as the system has acquired Algol characteristics with < Mg, 
the orbit widens and a liberal era can restart as soon as the mass- 
transfer-rate is sufficiently large. It happens that binaries run into 
superficial contact when during evolution M^ ^ Mg so that the 
impact parameter vanishes. This leads to a brief pause in spin-up 
and further construction of the hot spot on the gainer. After some 
time, when Mg exceeds the contact is broken, as suggested 
for the short period binary CN And by Van Hamme et al. ( 2001 ). 

Table 1 covers a wide range of initial periods for (9+5.4) M© 
binary at birth. Mass-loss from the system always occurs soon 
after RLOF ignition. At the beginning of this violent epoch, the 
binary does not meet the criterion of Peters (120011) stating that 
in a semi-detached binary, the donor needs to be less massive, 
cooler, fainter and larger than the gainer before the binary can 
be labeled as Algol. At the end of the liberal era the binary starts 
its life as an Algol. It is straightforward to determine the times 
that the Algol spends with changing values of mass-ratio and 
orbital period. 

From Table 1 we learn that short initial periods yield con- 
servative evolution only. The quantity of mass lost from the sys- 
tem rises with increasing initial orbital period. Binaries that start 
their RLOF at the end of hydrogen core burning (case- A) or at 
the beginning of hydrogen shell burning of the donor (case-B) 
lose a large amount of mass without being an Algol system for 
a long time. After a long era of Algolism during hydrogen core 
burning of the donor, RLOF and subsequent Algol status are also 
achieved during hydrogen shell burning. During this epoch the 
mass-transfer-rate is, however, never sufficiently large to trigger 
mass-loss from the system. 



initial orbital period of 3.6 days loses more than 3 M© during 
its phase of rapid RLOF, shows its Algol characteristics for ^ 4 
million years and liberal evolution yields a A ^ a little above 0.1. 
Since we limited ourselves in this paper to binaries undergoing 
RLOF-A at first, we plan to examine binaries with larger initial 
orbital periods so that the fast phase of RLOF-B occurs towards 
a system having an accretion disc. 

For our small sample of binaries with an initial primary mass 
of 9 Mq and q = 0.6, we used the initial orbital period dis- 
tribution of Popova et al. (1982) and compared the obtained 
mass-ratio distribution with the result from conservative evolu- 
tion. It is shown that the ^-bin [0.4-0.6] which was poorly pop- 
ulated through conservative binary evolution (less than 30%) 
is now represented much better (^60%). It is thus clear that 
our liberal scenario meets the observed ^-distribution of Algols 
better than conservative evolution. It is, however, doubtful that 
detailed overall correspondence between theory and observa- 
tions will be obtained, because all the low-mass binaries evolve 
almost conservatively, despite the fact that they also develop 
some spin-up and not so bright hot spots. Conservative evo- 
lution is, however, not the general rule for binaries and the 
more massive among them can lose a significant amount of 
mass into space. In a near future we will complete our website 
(http : //www . vub . ac . be/astrof ys/) containing an atlas of 
conservative evolutionary calculations with the results from the 
liberal calculations as presented in this paper. Future observa- 
tions of mass-transfer-rates, energy contents, surface areas and 
temperatures of hot spots will certainly enable the researcher to 
refine the evaluation of the quantity K for which relation ( 28 ) is 
only a first attempt. 
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9. Conclusions 

Mass-loss into the interstellar medium is possible during the 
short stage of fast mass-transfer soon after RLOF ignition when 
the binary is not yet an Algol. Low-mass binaries hardly show 
sufficiently large mass-transfer-rates to reach the critical rate 
needed to overcome the binding energy of the system, as re- 
quired by relation (21 ) for a direct impact system or (22) if the 
direct impact is succeeded by the formation of a hot spot at the 
outer edge of a Keplerian accretion disc. Spin-up of the gainer 
and hot spots are frequently created but the joint energy of both 
mechanisms will overshoot the binding energy of the system 
mainly for intermediate-mass binaries. Therefore we have ex- 
amined the evolution of the binary (9+5.4) M© with various ini- 
tial orbital periods. We calculate mass-loss into the interstellar 
medium for binaries with an orbital period larger than 1.8 days. 
We have calculated the amount of mass lost with respectively 
strong and weak tidal interaction and found that the influence 
of the tides on the mass-loss is not significant. Figure [7] shows 
the evolution of the mass-ratio q = ^ with time since the start 
of RLOF. The shortest initial orbital period (P = 1.8 d) evolves 
almost conservatively and yields a long living Algol that lives 
a very long time with q > OA. Table 1 indicates that the binary 
with an initial orbital period of 2.5 days loses ^ 1 M© during 
its phase of rapid RLOF. The system shows Algol characteris- 
tics for more than 10 million years and liberal evolution yields 
^-values which are systemically larger than those obtained with 
conservative evolution (A ^ a little below 0.1). A binary with an 
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Fig.6.yS as a function of time in and around the epoch of liberal (J3 < I) evolution, for a system of (9+5.4) M©, Pq = 3.2 d, for the 
weak and strong tide approximation. 





P=1.8 


F=2 


P=2.25 


P=2.5 


P=2.75 


P=3 


P=3.2 


P=3.6 


P=4 


strong tides 


AM 


0.236 


0.514 


0.889 


1.260 


1.811 


2.146 


2.540 


3.650 


4.011 


A tiiheral 


22,300 


20,900 


29,400 


31,500 


32,500 


29,700 


38,000 


32,200 


31,300 


A tAlgol 


9,858,800 


12,169,600 


13,817,500 


11,004,500 


8,949,000 


6,946,200 


5,746,800 


3,736,700 


27,800 


weak tides 


AM 


0.084 


0.023 


0.152 


0.933 


1.545 


2.391 


2.737 


3.315 


3.237 


^ liberal 


8,100 


10,500 


9,600 


11,700 


13,600 


19,300 


17,100 


18,700 


14,700 


A tAlgol 


10,067,000 


11,887,600 


14,570,400 


10,970,200 


8,778,700 


7,091,300 


5,970,500 


3,670,700 


1,412,600 



Table 1. Characteristics of liberal evolution of a binary (9+5.4) Mq undergoing RLOF-A. Systems with periods below 1.8 days 
evolve almost conservatively and are not mentioned. The mass-loss from the system hardly depends on the strength of the tidal 
interaction. The mass-loss is expressed in solar masses. The duration of liberal evolution and of the Algol phase is in years. 




Fig. 7. Evolution of the mass-ratio q = jf with time after onset of RLOF. Binaries on the thin part of a curve are no Algols, unlike 

those located on the thick parts. Curves are labeled as follows: (1) Conservative evolution with Pq = 3.6 d, (2) Liberal evolution 
with Po = 3.6 d, (3) Conservative evolution with with Pq = 2.5 d, (4) Liberal evolution with with Pq = 2.5 d, (5) The evolution with 
Po = 1.6 d which is always conservative. 
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